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Neutrino emission via the plasma process using the vertex 
formalism for QED in a strongly magnetized plasma is consid- 
ered. A new vertex function is introduced to include the axial 
vector part of the weak interaction. Our results are compared 
with previous calculations, and the effect of the axial vector 
coupling on neutrino emission is discussed. The contribution 
from the axial vector coupling can be of the same order as or 
greater than the vector vector coupling under certain plasma 
conditions. 
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I. INTRODUCTION 

The decay of electromagnetic oscillations in a plasma 
into neutrinos is of interest as a stellar energy loss mecha- 
nism |ffl,B . The presence of a plasma allows the refractive 
index to be less than unity, which is necessary for the de- 
cay of a photon into a vv pair. The resulting "plasma 
process" for neutrino emission, has been studied in the 
presence of an unmagnetized plasma There were 

also calculations in which the background medium was 
taken to be the magnetized vacuum; however the refrac- 
tive index for the magnetized vacuum is always greater 
than unity, and it was assumed that the presence of some 
low density plasma could lead to an appropriate open- 
ing of phase space to allow the process to proceed [|| |) . 
Here we perform a consistent calculation, in which the 
background plasma is included explicitly and thus the 
kinematic condition is not an ad hoc addition. 

The decay depends on the properties of the waves, and 
a magnetized plasma can support a variety of natural 
wave modes. Canuto, Chiuderi and Chou |jl(],[ll) con- 
sidered the plasma process in a magnetized plasma and 
they considered several possible wave modes. However 
their analysis neglected the axial vector aspect of the 
weak interaction and they also did not use the exact 
electron wavefunctions in a magnetic field. These defi- 
ciencies raise doubts about the validity of their results at 
high magnetic field strengths. 

In this paper we calculate the amplitude for the decay 
7 — ► vv in a magnetized gas of electrons to 0{Gp). We 
avoid the weaknesses in p-|^, |Io| , |rT| , by using the exact 
electron wavefunctions in a magnetic field and including 
the effects of the background plasma. In Sec. II, the for- 
malism required to treat V-A interactions in a strongly 
magnetized plasma is summarized and extended. The 
formalism used is the vertex formalism [|l^-|T4[], which 
allows both a momentum space representation for QED 



in a strong magnetic field, and a means to calculate the 
response tensors of a magnetized plasma. In Sec. Ill the 
transition rate for the decay of a given wave mode is cal- 
culated. Wave modes for a plasma with a cold electron 
distribution and for a thermal electron distribution are 
considered. It is shown that the results of may 
be recovered with suitable approximations. The implica- 
tions of the axial vector aspect of the weak interaction 
are discussed, and the neutrino emission rates from dif- 
ferent plasma modes are compared. Natural units with 
h = c = 1 and Boltzmann's constant, k = 1, are used 
throughout, and only standard neutrinos are considered. 

We find that the presence of a magnetic field has very 
little effect on neutrino emission relative to that in an 
unmagnetized plasma except for magnetic fields close 
to the "critical" magnetic field strength B c = m 2 /e = 
4.41 x 10 9 T. At high magnetic field strengths there is 
an enhancement of neutrino emission due to a large pro- 
portion of electrons being present in their lowest Landau 
orbital. We derive a criterion for determining when the 
axial vector contribution is likely to be important for 
neutrino emission. 



II. V - A INTERACTIONS IN A MAGNETIZED 
PLASMA 

A. Vertex Formalism 

A systematic development of QED in a strong mag- 
netic field was presented by Melrose and Parle |l2|-[l4| . A 
summary of the electron wavefunctions and vertex func- 
tions is contained in Appendix A. The electron energy 
levels for a static background magnetic field of magni- 
tude B parallel to the 3-axis are 

£ q = (m 2 +p\ + 2neB)i, 

where {q} labels the set of quantum numbers which in- 
cludes the parallel momentum, and the Landau levels, 

n = 0, 1, 2, The ground state, n = 0, is a singlet state 

and the states n > are doubly degenerate due to two 
spin states. In the Landau Gauge, the vector potential 
is A(x) = (0, Bx,0). The electron wave functions ipq{x) 
are eigenfunctions of a spin operator and the magnetic 
moment operator is chosen, cf. [ |L5| : 

p, = ma — i7 1 erx (p + eA(x)), (1) 

where er denotes the Pauli spin matrices. This spin oper- 
ator commutes with both the Hamiltonian and radiative 
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corrections to the Hamiltonian, and its eigenfunctions 
have symmetry between electron and positron states. 
Here electron and positron states are labelled by e, elec- 
trons correspond to e = 1 and positrons to e = — 1. 

To allow a momentum representation, wc use a vertex 
function [7^(k)] M , defined in the following way fll^] : 



1 

v 



dxexpHk-x)^,(x) 7 ^(x), (2) 



where V is the normalization volume. The incom- 
ing electron has quantum numbers {e, q}, the outgo- 
ing electron has quantum numbers {e',q'} and the 
outgoing photon has 4- momentum — (w,k) = 
(lu, k± cos - !/;, k± sin ip, k»). From the definition of the ver- 
tex function, the following symmetry property is clear: 



(b^( k )]T = [7^(-k)r 



(3) 



Since we also wish to treat the axial vector component 
of the weak interaction, we define an axial vector (AV) 
vertex function similarly to Eq. (||), replacing 7 M with 
7 M 7 5 . The AV vertex function is denoted by [7^4 (k)] 5 . 

An electron - photon vertex corresponds to the stan- 
dard vertex function, and an electron - Z boson ver- 
tex corresponds to a combination of the standard vertex 
function and the AV vertex function. The V-A theory 
of weak interactions ignores Z and W boson propagators 
and considers only the charged and neutral currents at 
a point interaction. The neutral current component of 
the interaction may be expressed using the standard and 
AV vertex functions, and hence the charged current com- 
ponent of the interaction may also be expressed in this 
manner through the use of a Fierz transformation (see 
e.g. ). 

Another vertex function, [r^(k)] M , which is a gauge 

invariant part of [7^(k)] M , is identified due to the de- 
sirability of having a gauge invariant theory. The gauge 
invariant part of the AV vertex function is identified as 
[r^(k)]g. (An example of the separation of gauge de- 
pendent and gauge independent terms is given in Ap- 
pendix A). 

Using the vertex formalism, we obtain a momentum 
space representation of the effective V-A interaction La- 
grangian as: 



G p 1 c: 

Arff = --/=u(gi)7 M (l - 7 M<?2) 



{^l[r^(k)]^ + B[r^(k)]^}, 



(4) 



where e, q and e', q' label the incoming and outgoing elec- 
tron states, respectively, and u(qi) and v{q2) are the neu- 
trino and antineutrino wavefunctions, respectively. The 
Fermi constant is represented by Gf and the constants 
A and B are given by A = 2 sin 2 9w + \ and B — — \ for 
electron neutrinos, and by A — 2 sin 2 6\y — \ and B = \ 
for muon and tau neutrinos, where #w is the Weinberg 



angle. In the approximation where sin 2 9w = \, then 
.4 = for muon and tau neutrinos, so that only the axial 
vector component of the weak interaction contributes to 
their emission, as noted by 

The gauge invariant form of the AV vertex function 
is given in Appendix A. When calculated using magnetic 
moment operator eigenfunctions, it obeys similar symme- 
try relations to the standard vertex function. Explicitly 
these relations are 



([r£(k)]2r = [r£(-k)]S, 

[l7,'- £ (-k)]£ = (-l)''-^4(k)]£. 



(5) 
(6) 



B. Response Tensors 

The standard linear response tensor for a plasma, 
II A " y (fc) when written in covariant notation satisfies the 
equation 



(7) 



where j M (fc) is the induced 4-current and A fi (k) is the 
fluctuating part of the electromagnetic field. In the ab- 
sence of a plasma, II Atl/ (/c) is the vacuum polarization 
tensor. Using the vertex formalism, one can introduce 
a medium using the electron occupation numbers, and 
then with the assumption that the occupation numbers 
of a state are independent of spin, the linear response 
tensor in a magnetic field becomes ( [fill ): 



n',n—0 e',e— ± 



^||{|(e'-e) + e^-e'^} Tr! 



2tt uj — e£ q + e'E q ' + iO 



where Ttf e is the product of vertex functions summed 



over spin states, i.e. 



e'e 



flu 



a' ,er— ± 



[r-(k)nr"(k)f 



(9) 



The result of the summation in Eq. (^|) is presented in 
Appendix A, the 3-tensor form of the response tensor 
was calculated in |l7|] , correcting the result of , and 
the renormalized vacuum polarization tensor has been 
treated using the vertex formalism The infinitesi- 

mal imaginary term in the denominator of Eq. (^) arises 
from the requirement that the response tensor be a causal 
function. Note also that conservation of momentum is 
implicit through the relation e'p'^ — ej»|| — fc|| . 

The matrix element for the decay of a photon into a 
neutrino pair contains the product of a standard vertex 
function and an AV vertex function. This allows one 
to identify an axial vector response function B, which 
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can be generalized in the same manner as the vacuum 
polarization tensor to include a medium. Thus we have 



2ir ^ ^ 

n' ,7i—0 e', e— ± 



2vr w - e£ 9 + e'£ q , + iO 5 e ' e ' 1 j 
with a similar sum over spin states 

bT £ = [r ^(k)]^[r^(k)]r. (ii) 

cr' , CT— i 

The axial vector response for a magnetized plasma 
has not been written down previously to the best of our 
knowledge, however there have been calculations of tt 11 ^ 
using the proper time formalism for the magnetized vac- 
uum m ■ The result of the sum over spin states is pre- 
sented in Appendix A. 

III. MATRIX ELEMENT 

There are two contributing diagrams to the amplitude 
for the plasma process for neutrino emission to 0(G 2 F ). 
These diagrams, and the diagram for the process when 
regarded as a V-A interaction are shown in Fig. 1. For 
simplicity, only electron neutrinos are considered here. 
Using the V-A theory of the weak interaction, the matrix 
elements for the W and Z diagrams contributing to the 
decay may be expressed in the same form using a Fierz 
transformation. The matrix clement is 



M 



G, 



fi 



-u(qi)-f"(l - ~f 5 )v(q 2 ) 



V2e 

x (AU^(k)+BUjl(k))A^k), 



(12) 



where is the fluctuating part of the electromagnetic 
field. For a magnetized plasma, the only difference in the 
matrix element Eq. ( |l2| ) from that for an unmagnetized 
plasma is the form of the electron wavefunctions. 

The decay rate for a plasma mode is taken to be the 
transition probability per unit volume of x — k space 
per unit time for the decay of a quantum of a plasma 
mode into a neutrino antineutrino pair. Provided that 
the refractive index of the plasma is close to unity then 
the decay rate T m for a mode M can be written in the 
well known form (e.g. |,|,@): 



M 



G 2 F R M (k) 
&ire 2 e \u)M{k)\ 



(13) 



where i?M(fc) is the ratio of the electric energy to total 
energy in the wave, and 



Q^ = \{An <Til {k)+BIl <T l{k))el I {k)\* 
x \(An TV (k)+ BU T l{k))e T M {k)\. 



The power emitted per unit volume for a given wave 
mode may be obtained by multiplying the decay rate for a 
given mode by the energy, ujm, of a photon in that mode 
and the plasmon occupation number /(ujm)- Integrating 
over momentum space leads to the power emitted per 
unit volume by the decay of wave quanta into neutrinos 
as 



d 3 k 

(2tt) 3 



wmFm/(wm)- 



(15) 



To calculate the neutrino emission from given plasma 
conditions, one proceeds through the following steps. 
First, the electron distribution function is required to 
determine the plasma response. Secondly, the plasma 
response is used to find the natural wave modes of the 
plasma. The combined effects on the wave properties of 
the vacuum polarization of the magnetized vacuum and 
of the plasma response was discussed by [plf ; but the 
vacuum contribution is ignored here. Thirdly, the po- 
larization vector, dispersion relation and response tensor 
are used to calculate the decay rate, Eq. (13), which is 
then integrated in Eq. ( |l5| ) to determine the power emit- 
ted in neutrinos. The most difficult step in obtaining 
analytic results is solving for the wave modes in a magne- 
tized plasma - there are relatively few cases in which the 
modes are simple enough to allow computational ease. 
However, given the relatively straightforward procedure 
to calculate the power emitted in neutrinos, there is the 
opportunity to obtain numerical rates for a large range 
of plasma conditions. 

A further simplification for computational ease is to 
take the long wavelength limit of the expressions for W LV 
and Il^g. These expressions are presented in the Ap- 
pendix. Rather than using the response 4-tensor to de- 
termine the wave modes, it is convenient to use the di- 
electric 3-tensor which is related to the response 3-tensor 
by 



1 



rlT 



Ki j = 5i j+ ■ , 

The form of the dielectric tensor is the same as that for 
a cold plasma [E2| 



K\ 



S -iD 
iD S 
OOP 



(16) 



A. Neutrino Emission from a "cold" plasma 

Starting from Eq. ( |l6| ) for the dielectric tensor for a 
cold plasma, the equation for the refractive index n takes 
the form 



(14) 



An A - Bn 2 + C = 0, 
with the coefficients 



(17) 



3 



A = S sin 2 9 + P cos 2 9, 

B = (S 2 - D 2 ) sin 2 9 + SP(1 + cos 2 6), 

C = P(S 2 — D 2 ). (18) 



A specific solution n — um of Eq. (|17]) defines the mode 
M. The general expression for the polarization vector of 
a mode may be expressed as p0[]: 



Kmk + T M t + ia 

(K 2 M + T 2 M + l)^ 



(19) 



where k, t and a are given by 



k = (sin#, 0, cos#), t = (cos 9, 0, — sin0), 

a =(0,1,0), (20) 

and the coefficents Km and Tm are given by 



K M 



'■M 



)£>sin( 



An 2 M 



PS 



DP cos 9 



M 



Mi 



PS 



To calculate the natural modes from the dielectric 
tensor, the wavevector is taken to be k = |k|/c = 
|k| (sin 9, 0, cos#), where 9 is the angle between the 
wavevector and the magnetic field. Note that the choice 
of gauge here is the temporal gauge, thus the polarization 
vector for a mode M takes the form e M = (0, e M ), where 
en is the polarization 3- vector. The polarization vectors 
take simple forms for the cases 9 = and 9 = ir/2. For 
9 = 0, there can be two circularly polarized modes (or 
only one if the other is evanescent), the ordinary and 
extraordinary modes, or only one longitudinal mode, the 
others being evanescent. For more general angles of prop- 
agation, i.e. 9 7^ 0, 7r/2 the modes have neither purely 
longitudinal or purely transverse polarization. The dis- 
persion relations and polarization vectors for the modes 
at 9 = and 9 = ir/2 are given in Table 1. 

A cold plasma electron distribution is 



fn i^n 



4^ 
~eB 



n + +n ] 5 n0 S(p\\), 



(21) 



where the n £ (£„) correspond to the number densities of 
electrons and positrons. Equation (|2l]) is a distribution 
in which all the electrons are in their lowest Landau or- 
bital. Having substituted Eq. (|2l]) into the expression 
for the dielectric tensor, if we then take the classical limit 
(lo <C m, eB <C m 2 ) and assume a purely electron plasma, 
we obtain the dielectric tensor of magnetoionic theory, 
which is exactly that used in previous investigations of 
the plasma process in a magnetized plasma [0,0 . 

Using the notation X = u> 2 /oj 2 , Y = Q e /ui, where ui p is 
the plasma frequency (defined by o> 2 = e 2 n e /eom ) and 
£! e is the electron cyclotron frequency, the magnetoionic 
theory gives GO): 



X 



iXY 



-Y 2 

K 3 3 = l- X, K l } = otherwise. 



, (22) 



1. Emission at 9 = 

The longitudinal mode at 9 = is independent of the 
magnetic field and the decay rate is the same as the 
known value for an unmagnetized plasma [^| ^) . The 
power emitted in the transverse modes at 9 = can 
be written in a particularly simple way. Taking A = 1 
to label the ordinary mode and A = — 1 to label the ex- 
traordinary mode, the power emitted per unit solid angle 
is 



Q> 



G 2 F 



3847r 5 a 



^Ajfl-n^/H (23) 



where n\ is the refractive index for the mode A. The 
frequencies Lu min and w max correspond to the frequencies 
at which the refractive index is and 1 respectively - for 
a cold classical plasma w max = oo. The frequencies here 
satisfy ui -€1 m, so this value is used for the upper cutoff. 
Integrating from Lu m - ln might appear to contradict the as- 
sumption made in deriving Eq. (p^), that the refractive 
index be close to unity. However, for almost all of the 
frequency range the refractive index is close to unity and 
thus the results obtained here are not compromised by 
this — it is a far less serious approximation than the as- 
sumption of a cold plasma. For a plasma described by 
the magnetoionic theory, Eq. (^3|) reproduces the results 

of®- 



2. Neutrino Emission from modes at = tv/2 

The behavior of the ordinary and extraordinary modes 
at 9 = ir/2 is less simple than for 9 = 0. The axial vector 
part of the weak interaction can only couple to modes 
which have a component of their polarization vector par- 
allel to the magnetic field. The ordinary mode has such 
a component but the extraordinary mode does not. For 
comparison, the power emitted by the ordinary mode per 
unit volume per unit solid angle is presented both with 
and without the axial vector coupling. With the axial 
vector coupling one has 



Gp 



15367r 5 a 



du w s n (l-<) 2 (4-3n 2 )/( W ) 



(24) 



and without the axial vector coupling one has 



Qo 



G 2 F 



dto w 8 n„(l - njr/H- (25) 



384?r 5 a 

The power emitted in the extraordinary mode is 

du) uj s n x (\ — n 2 ) 
4D 2 S(S-1 



Gp 



384?r 5 a 

(S'-l) 2 + L» 2 



D 2 + S 2 



(26) 

/H- (27) 
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As can be seen in Fig. 2, the inclusion of the axial vec- 
tor coupling leads to results which differ significantly 
from those obtained when it is ignored. The effects are 
most pronounced for lower electron number densities, 
and there is only about a 25% increase in emission close 
to the peak when AV effects are included. 

Using the dispersion relation Eq. ( fL7| ) the plasma res- 
onances (which correspond to n 2 — oo) may be found 
from the equation A/C = 0. In magnetoionic theory, the 
solutions to this equation are 

-1(0) = \H + Hi) ± \{H + - H^l cos 2 6}i. 

(28) 

Canuto et al jflj claimed that w 2 (0) = £1 2 is a mode 
which can lead to enhanced neutrino emission at high 
plasma densities (p > 10 11 g cm -3 ). However, the plasma 
resonance does not satisfy the kinematic condition that 
the refractive index be less than unity, which is required 
for the plasma process to proceed — hence no energy 
can be lost through this mechanism. We conclude that 
there is no such enhanced emission at exceptionally high 
plasma densities. 



B. Neutrino Emission from a Thermal plasma 

To obtain analytic expressions for the energy loss in 
neutrinos from a thermal magnetized plasma one assumes 
a thermal form for the electron distribution in the expres- 
sion for the response tensor Eq. (||). We make either a 
non-relativistic or semi-relativistic expansion of the reso- 
nant denominator in the tensor to simplify the analysis. 



1. Role of the axial vector coupling 

For an unmagnetized plasma, it has been shown nu- 
merically that the AV contribution to energy loss via the 
plasma process is of the order of 0.01% for temperatures 
below 10 11 K | p3| . However, in a magnetized plasma, it is 
possible that the AV coupling can have a more significant 
effect on neutrino emission. Physically this may be seen 
as follows: the AV coupling cannot affect processes in a 
system which has reflection symmetry; it requires that 
there be some axial vector in the system to which it can 
couple. Although there is no such axial vector in a clas- 
sical magnetized plasma (apart from the magnetic field, 
which is not relevant here), in a quantum magnetized 
plasma, the electronic ground state (the lowest Landau 
orbital) corresponds to a specific spin state, unlike all 
excited states which have two degenerate spin states. A 
plasma with a significant fraction of its electrons in their 
lowest Landau orbitals thus has an appropriate axial vec- 
tor that allows coupling to occur. 

Hence we expect the AV component of the weak in- 
teraction to be important when a significant fraction of 



the electrons are in their lowest Landau level. Consider 
a Fermi distribution of electrons 



9n 



exp[(£ q -p)/T\ + l' 



(29) 



where g n is the degeneracy of the nth energy level, /i is 
the chemical potential and T is the temperature. Tak- 
ing the limit that T becomes large, Eq. (|2^) becomes a 
Maxwell - Boltzmann distribution, 

f(£ q ) = g n exp[-(£ q - n)/T], 

which when normalized in the non-relativistic limit gives 



f(£ q 



■ 9n- 



An 2 n e tanh(A/2) 
eB (2tttoT)2 



exp 



2toT 



nX 



(30) 



where A = eB/mT. The normalization that is used is 



oo „ 

/ d P\\9nf(£q) = 
n=0 ^ 



47r 2 n e 
eB 



(31) 



The most probable value of n should be when the en- 
ergy associated with parallel motion is of the same or- 
der as that associated with perpendicular motion, corre- 
sponding to equipartition of energy. This occurs when 
p 2 /2mT = nX. The condition is 



F 



neB 



(32) 



2mT mT 

If one replaces p| by (p 2 ) — mT and n by (n), Eq. ( |3^ ) 
gives 



muT 



2ehB 2ekB 



0.06§, 



(33) 



where T is the temperature in kelvin and B is the mag- 
netic field in tcsla. For a young, highly magnetized 
white dwarf star, one might expect the surface values 
of magnetic field and temperature to be B ~ 10 5 T and 
T < 10 6 K respectively, which gives (n) ~ 0.6, so a size- 
able proportion of the electrons are in their lowest Lan- 
dau orbital. Hot, strongly magnetized white dwarfs or 
their precursors, and neutron stars, are objects which 
are likely to have their plasma neutrino emission affected 
by the presence of a strong magnetic field. 

If we take (n) < 1 to characterize when most of the 
electrons are in their lowest Landau orbital, then a cri- 
terion for whether the AV part of the weak interaction is 
important for neutrino emission is 



> 0.06 



(34) 



Temperature and magnetic field regimes that occur in 
the interior of neutron stars and white dwarfs are com- 
pared with the criterion Eq. Q) in Fig. 3. Whilst there 
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is a larger range of B and T conditions available for neu- 
tron stars, the electrons are almost certainly degenerate, 
in which case, the refractive index is greater than unity, 
and the plasma process is forbidden. The white dwarf 
and neutron star internal conditions are taken from |jj . 



and for cyclotron absorption we have 



to — £ q + £ q i ~ (u> — Q e ) 



P\\ k \\ 



m 



The resonant response tensor becomes 



(39) 



2. Neutrino Emission 

To calculate the response tensor, start with the distri- 
bution function Eq. ( p(i| ) and substitute into Eq. (||). In 
the absence of positrons, the response tensor takes the 
form 



n""(fc) = 




■ £ q + £ q > + iO ++ 



(35) 



may be found by replacing by sTfi . 
In making the non-relativistic approximation, it is as- 
sumed that the thermal energy of the electrons is much 
less than their rest mass energy, i.e. T « ra ~ 6 x 10 9 K 
plf . Provided that one considers a plasma with VL e <C m 
and modes such that to < Q e , so that £ q + £ q i > 2m 3> u>, 
the second and third denominators do not vanish. There 
are three transitions that can occur in the highly mag- 
netised plasma, all of which must be taken into account 
when calculating the plasma response. There are pro- 
cesses in which an electron remains in the same Landau 
orbital after emission, i.e. n = 0, n' = 0, and there is 
also cyclotron emission (n = l,n' = 0) and cyclotron 
absorption (n = 0, n' = 1). The first denominator has 
a zero, called a resonance, corresponding to either cy- 
clotron emission or cyclotron absorption. The resonant 
term is sensitive to finite temperature effects, but the two 
non-resonant terms are not [|24| . Hence for the two non- 
resonant terms, one may set the distribution function to 
be 



f(£ q ) 



47r 2 n e 
eB 



S(p\\)SnO, 



(36) 



These "non-resonant" contributions to the response ten- 
sors are shown in Appendix B. 

In the resonant terms we use a Maxwellian distribution 
of electrons in their lowest Landau orbital: 



i 



47r 2 n f 



y/2nmT eB 



Sno exp 



2mT 



(37) 



For cyclotron emission the resonant denominator is 

lo - £ q + Eg, ~ (w + fi e ) - 



(38) 



n^(fc) 



dpn exp 



n 



r^(n = 0,n' = l) T^(n = l,n' 



0) 



(u> + O e ) — pukn/m {uj — Q e ) — P\\k\\/i 



with the expression for n AI 5 rcs obtained by replacing 
by hT^. This leads naturally to expressing the compo- 
nents of the tensors in terms of the plasma dispersion 
function 4>{z), defined by pof 



dt 



t 



(42) 



where z = (u> ± fl e )^/m/k\\ V 2T. To simplify the expres- 
sions obtained in Appendix B further, one can make use 
of the asymptotic expansion of 4>(z) for large z EQ] 



1 



1 

2^2 



3 
Az 1 



'irze 



(43) 



taking only the highest order term in z and ignoring the 



imaginary part. This leads to Eq. (B10) for the response 
tensor. For waves propagating parallel to the magnetic 
field, the tensor reduces to the same form as for a cold 
plasma; the power emitted in the ordinary and extraor- 
dinary modes is given by Eq. (p3|), and the longitudinal 
mode is identical to that for an unmagnetized plasma. 
Waves which are not propagating parallel to the mag- 
netic field, lead to expressions for power emission which 
are more complicated than Eqs. (]23|)-(|25|) and may have 
significant contributions from the axial vector coupling. 



IV. CONCLUSIONS 

The calculations in this paper address a number of is- 
sues relating to the neutrino plasma process. Firstly, the 
plasma process for neutrino emission has not previously 
been calculated taking into account strong magnetic field 
effects, plasma effects and the axial vector part of the 
weak interaction. Our work takes a consistent approach 
to the inclusion of a plasma and the kinematic conditions 
under which the neutrino plasma emission process may 
occur, as opposed to the inconsistencies in previous calcu- 
lations • Secondly, the formalism for looking at weak 
processes in a strongly magnetized plasma, has not previ- 
ously been able to deal with diagrams containing electron 
loops. The axial vector vertex function and axial vec- 
tor response tensor described here provide mathematical 
tools which can be used for such calculations. Thirdly 



G 



we have produced some analytic approximations to the 
power emitted in neutrinos from a volume of plasma with 
a given magnetic field, electron number density and tem- 
perature. The exact results for the response tensors mean 
that these can be used to calculate numerical results for 
magnetic fields greater than the critical magnetic field 
- this regime has not been investigated here. We have 
also shown that contrary to the case of an unmagnetised 
plasma, the axial vector coupling can have a role in af- 
fecting neutrino emission via the plasma process, and we 
have suggested a simple criterion with which to estimate 
whether such axial vector effects are likely to be impor- 
tant. 

The neutrino plasma process is related to neutrino 
Cerenkov radiation by a crossing symmetry (see ||), so 
that the results obtained here for the response functions 
can be used to study the Cerenkov process in plasmas 
with a refractive index greater than unity. 

The magnetic field dependence of the plasma process 
parallels the results found for neutrino dispersion in a 
strong field p5[ , in that the results are relatively insen- 
sitive to the magnetic field. The plasma process is only 
sensitive to the magnetic field for B close to B c - there is 
a much stronger dependence on temperature and electron 
number density than the magnetic field. Hence unless 
one considers strongly magnetized plasmas, most of the 
expressions derived for unmagnetized plasmas are ade- 
quate. 

There are several highly magnetized astrophysical en- 
vironments where the plasma process may be of impor- 
tance. These are in the cooling of giant stars with highly 
magnetized cores, in the early stages of the evolution of 
a hot magnetized white dwarf. The process may also be 
of importance for neutrino emission from neutron stars. 
The enhanced neutrino emission due to the axial vector 
coupling in regions of stronger than average magnetic 
field might contribute to an anisotropic neutrino lumi- 
nosity which has been suggested as a possible mechanism 
for the large proper motions of many pulsars. 
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APPENDIX A: 



The electron wavefunctions determined using the magnetic moment operator have been determined by |L2| |, we 
display them here for convenience 



^(x,t) 



-ieS Q t 



where 



e~xp(iep y y + iep\\z) 
{4£,£0(£,+£0)( f o +m)} * 



->e,-l 



5a,l 



S<j,l 



/(£, + £°)(£°+m)v,,_i(0\ 

-iPnP\\V n (t) 

P||(f° + m)v n _i(0 

\ ^„(^+f 3 °) V „(e) y 

P||(£S+mK_i(0 
-ip„(f, + fijK(0 
(f, + £j)(£j+mK_i(0 
*PnP||Vn(£) 



+ 5 a ,-l 



( 



\ 



+ s a 



/ -iPnP\\Vn-l(0 \ 

~ip n {£ q 
\ - P ||(£ 9 + m) Wn (0 / 

/ ip„(£,+£ g °)u Il _ 1 (o y 

-P|l(££ + mK(0 
ip n P||f„_i(0 
V(£,+£ g )(£,°+™K(0/. 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 



where the quantities p n ,£ q and £ q are given by: 



(2neB)-*, £° q = (m 2 + p 2 n Y, £ q = {£^+pp 



(A6) 



In Eq. (A.2), er is the spin quantum number which takes the values ±1 for spin up and spin down respectively, and e 
is the sign of the energy. If p^ is the z component of momentum for an electron (e = 1), then pii represents minus the 
z component of momentum for a positron (e = —1)- The functions v n (£) are normalised simple harmonic oscillator 
wavefunctions of the form 



MO 



tf n (Oexp(-k 2 ) 



( 7r l/2 2 n n! )l/2 ' 

where H n (^) is the nth Hermite polynomial and 

Z = (eBy/ 2 (x + ep y /eB). 

The separation of the vertex function into gauge dependent and gauge independent terms was given in ]T^ | for 
several choices of electromagnetic gauge. We write down their result for the Landau gauge 



[Y q , q (W = {(2Tr) 2 /V(eB)i}exp[ik x (ePy + e'p' y )/2eB]S(ep y - e'p' y - k y )S(ep z - e'p' z - fc z )[r^(k)]^. 



(A7) 



A similar separation may be made for the AV vertex function. One may write out the gauge invariant vertex function 
for the magnetic moment operator eigenfunctions as: 

[r^(k)f = (%C g [8M<4\(4-l + Pn'PnJl^l 

e^' e ? (-Pnexp(iCT?/')4t-1_ CT - P' n > exp(-iai))jj,_ l+a ), 
ieaf3 qlq (p n exp(iaip)J l +_ a l _ lJ - p' n , exp(-iaip)J l v _ l+a .), 

-e(j5^- a {a e q , e q {-p n e^p{ia^)J l +_ a l _ ty + p' n , exp(icr^) J;',_;_ CT ), 



%eab^\{j\,-l + P'n>PneM2*^)J l /-i-2*)> 
d q > £ q (-Pn exp(iat/)) 4t^_ CT + p' n > exp(icrV0 J i'-;- CT )}]> 



(A8) 
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where the argument of the J functions is the same as in Eq. (A14), and with the coefficients: 



and the abbreviations: 



a q'q 


= 


-p\\p\\) 


+ 0-5 e '_ e (P|| 


+ P|) 


Pq'q 


= Wi- 


-p\\p\\) 


4- cr<5 e /„ e (p'| 


-pi) 


4\ 


= 


+ p\\) + 






e'e 

a q >q 




+ p\\)- 


a5 e >- e (l + 


p||P|) 


/ 

°q>q 




-p\\)~ 


(T^_ e (l- 


0||P|) 


a q'q 




-p\\p\\) 


- <7<5 £ '-e(p{| 


+ p\) 



(A9) 



C 9 = 



'(5 9 + £: 9 )(£:° + m)^ 



^£q£q 



' (£ q ,+£ q ,)(£° q ,+m) 



{iexp(iip)} 1 , 

1 

2 

{iexp(iV0} J > 



(AfO) 



P|| = Pll /(£q+£° q ), 

p\ l =p\ l /(£ q/ +£ q ,), 

Pn =Pn/(£° +m), 

p' nl =p n -/{£° ql +m). (Aff) 

As noted by |T^ ], the gauge invariant vertex function satisfies symmetry properties similar to Eq. (^|), provided that 
the magnetic moment operator eigenfunctions are used. These symmetry properties are: 

([r^(k)f)* = [r-:(-k)f, (Ai2) 
[r-,;'- e (-kr = R''- ; [r^(k)r (ais) 



with the symmetry property Eq. ( Af 3 ) holding for suitably chosen phase factors. 
The gauge invariant form of the AV vertex function is 



[r^M = Cg^W^Ml'-i - Pn'PnAtV, 

Kq,q{-Pn eX P (icn/')4-1- ( T + Pn> eX P(-* Cr V')4-i+ C r)> 

ieair qlq {p n exp(iaip) J l ^"i_ a + p' n , exp(-io-^)j/,_ i+CT ), 

O q >q(J l l>-l - Pn' PnJ l V-l)} 

-eaSv'-vif^qi-pnexpiia^J 1 ^^ - p' n , exp(iatp) J l V -i- a ), 



tf/M'-l + Pn'PneM2i^)J l +_ a l _ 2a ), 
iz09i'q{J l V-l ~ Pn'Pn eSXp{2i(T1p) J l /_^_ 2g ) , 

hp q (- Pn exp(i(rV)4tl_ CT - Pn' exp(i<"l>)4-i-o)}], ( Al4 ) 

where I = n — i(cr + 1) and a = ±1 is the spin eigenvalue (note that for the ground state spin singlet n = 0, the 
spin eigenvalue is a = —1). The J functions have argument k±/2eB, and are related to the generalized Laguerre 
polynomials (e.g. (29]) via the relation 

Mx) = ( (^yy) 2 e^/V/^Or) = (-)" JI!+». (A15) 







The properties of these functions have been summarized previously Jl2| . The coefficients for the AV vertex function 
are 
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<rS e 'e(p\ + P\\)+ <5e'-e(l + P\\P\\), 
(r$e'e(p\ ~P\\) + S e f- € (1 - p\p\\), 
CT<5 e 'e(l + + 5c'-e{p\\ + P\\), 

f}\ = -^.(l + P\\P\\) + Se'-e(p\\ +p\\), 
9q' q = -vSe'eQ- - P\\P\\) + 5 e '-e{p\\ ~ P\\), 

h q % = -*6 e , e (p\ + + ^_ e (l + p\ lP]] ). (A16) 



1. Response Tensors 

The sum over spin states for the linear response tensor and the vector axial response are presented below: 



T 



if (m»±pj |P ||) | [(j „ )2 + (j _ 1)2] ± ^ lj „ jrlj 

Z I Cq'Cq I Cq'Cq 

mil _ 1 f 1 ^ ( m2± P||P|l) l „ tt. ^3 ■ / Tn-1^21 , PnJV «! ™-l 

J ±££ — 9 S 1 T c? c? ( HVlJ + l J i/+l 7 J ± c c J v-l J u+l ) 

Z I Cq'Cq I Cq'Cq 



1 



T >22 ' J -, _( m "'-^'ll^'ll) 1 r/ T n \2 , / Tn-l\2l _ -PnPn' m jn-1 

1 ±ee — Ty \ 1 + J~~c ( W>v-\) + l-V+1 ) I + c c J v-l J u+l J 

Z I Cq'Cq I Cq'Cq 

Z I Cq'Cq I Cq'Cq 

rpOl J / tti-1 m I ni tti- 1\ _i_ jn jn , m- 1 xn — 1\ 

_ l jn jn _ jn—1 jrt— l 1 



rp\)Z J £_"■/ jn—1 jn _ jn jn—X\ i *——( T n 7 n _ 7 n— 1 7" n— 1 ^ 



i { (m 2 ±p\ m ) s \ 
Til = - J 1 T g T' K^i 1 ) 2 - (^-i) 2 )]- 



2 I £g 



13 _ _£ J ^"^11 r jn jn-1 , 771-1 mi , Pn'P\\ r ra pi , rn-1 TTl-ll 
tee — r, S c c J y "T" J i/ J 31 c c "+" J u J u+1 J 

Z CnCn' CnCa' 



rp23 _ _* e J r 7 n m-1 _ 7 n-l nil i Pn'P\\ r ™ rn tti-1 Tn-li 1 

J ±ee — r, S c c V J v-\' J v J u+1 J v\- S1 c c V J V J u-l J v J v+1 \ ( • 

I CqCq' CqCq' 

where v — n' — n. The remaining components may be constructed from the Onsager relations, which embody the 
requirements of time reversibility, and for a plasma with a static background magnetic field may be written in the 
form 

n°V,-k)|_ Bn =n°V,k)| Bo , n°>,-k)|_ Bo = -irV,k)| Bo , * 
ntf(w,-k)|_B = n*( w ,k)| BD . [ALii) 

The components of the 3-tensor part of may be constructed using the Onsager relations from the components 
given below, and the sum over spin states gives: 



10 



rpll 



5^± 



rp02 _ 
J 4-« — 



5^± 



rp03 



j.12 _ 



rril3 
5J- 



1-1- 



2:i 



T 22 



1 

2 I Sq 



rpUl _ ~ 
J +cc — T 



[(■C-i) -WTiT], 



r tti jn—1 _ Tn— 1 tiz ] -,- PjjjP^ r jn — 1 jn—1 jn j7i ] I 

l J v J i/+l J y v—l\ \ c c V J v J v+l v v — lj 1 ? 



2 I Sq 



f ,f 



ie \ P\\Pn _ x _ x P\\Pn' r n _x n _x n n A 

' \ F F ^ v v+ l v ^v—XX ^ c c V*v J v+1 J v ,J v-\\ ( 5 
Z I Cq'Oq Oq'Cq I 




i|^-g][( J "-i) 2 + (« 7 " + "i) 2 ]= 

£ J r 7« — 1 jn Tii m- li Pn' r 7-n— 1 Tn— 1 m m 1 

2[£ ^"+1 " J u-l J v Jt ~£~,\- v J "+l ■ / f« y K-:U 

Tn— 1 rri I m Tn — 1] -,- P n ' r Tn— 1 Tn— 1 i Tn Tn 
c- [■'/■- I. •',< '"' -V-rv It g^v 1 "V 'V l 



rri30 

5-t- 



2" 



-03 



Most of the components of the tensor Ir 5 satisfy the Onsager relations in the form of Eq. ( |A18 
3-tensor components and II g and II |. However, II | does not exhibit the symmetry in Eq. (IA1S 



(A19) 

specifically, the 
Physically, the 

reason for this is that the axial vector response violates parity. If the failure to satisfy the Onsager relations were due 
to the non-conservation of the axial vector current, one would expect that the components to be affected would be 
associated with k v . This is not the case, so we can ascribe the failure to satisfy Eq. flAl8f) to parity violation. The 

one also appeals to 



IS 



Onsager relations are derived from time reversibility, but to write them in the form Eq. 
arity and the reality condition for Fourier transforms. The tensor iE'j was calculated for a magnetized vacuum by 
and we note that when Eq.(22) of M is written in the co-ordinate system used here, it fails to satisfy the Onsager 
relations only for the II | and n 3< 5 components. 



2. Plasma response in the long wavelength limit 

When one takes the long wavelength limit of Eq. (||) , the linear response tensor becomes (cf. jij 

+ (f S | fS\ £n + l+£n ( U <> ma +*$> )) 

Un+1 In ] ^ - (£ n+1 + £ n f \ + S n+1 £ n ) f ' 

2vT 2^} U " +1 Jn > - (£ n+1 - £ n )2 [ £ n+1 £ n J 

(to 2 +p\)\ 1 



, (f D _ f D ] ^ -, 

Un+1 Jn) u 2 -{£ n+1 +£ n ) 2 \ ' £ n+ i£ 



2tt 2tt 1 Jn+i u; 2 -4£? l+1 \ £l +l ) Jn ^~A£l \ ft 



n=0 

n 00 = jjOl = n 02 = n 03 = n 13 = n 23 = ^ (A2Q) 

where f% = n + {£ n ) + n~(£ n ),f® = n + (£ n ) — n~(£ n ) and £ n = (m 2 + v\ + 2neS)^, noting that n + is the electron 
occupation number and nT is the positron occupation number. 



11 



Application of the long wavelength limit to yields 




(A21) 



By inspection, it is clear that if the electron and positron distribution functions are even functions of pii then the IT 1 ^ 
and II 2 5 terms are identically zero. Considering the other two components, n°| and n 3 °, rearrangement of the sum 
over n shows that only the n = term can contribute. The n = term vanishes identically for n°|, which means 
that the only non-zero component of the axial vector response tensor is n 3 ° . 



APPENDIX B: 



The results required to determine the response of a thermal magnetized plasma are summarized below. The only 
J functions that are required under the assumptions made in Section III are Jg and J®. These have the simple forms 



J °(w) = e-i u , 



(Bl) 
(B2) 



which are used below. The "non-resonant" contributions to the response tensor are 
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(B3) 



and for the axial vector response, 
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(B4) 



Due to the simplifications made in assuming a non-relativistic non-degenerate plasma, only one type of integral needs 
to be evaluated, namely 



h = 



dp\\ 



P 



a± - bp\ 



exp 



2mT 



I = 0,1 or 2, 



(B5) 



where T is the temperature, a± = uj ± f2 e and b — ku/m. The integral in Eq. (B5) may be re-expressed in terms of 

ich allows 

v27rmT 7/ N r \/2nmT , 



kj/r 



the plasma dispersion function <^(z) defined in Eq. (143), which allows one to obtain the following forms for the J; 



(z±), h = 



-M*t)-i), 



a±\/2nmT - 
h = ^ (<P(z±) - 1), 



(B6) 
(B7) 



where z± = a ±jb\ /2mT . 

Using Eq. (A17), the expressions in the resonant part of the response tensor may be evaluated 
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and for the axial vector response 
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Finally, we display the result for when one makes a high z± expansion, corresponding to J 1 « m. (The result 
for flg^ may be calculated similarly). 



n 00 = 



n 11 
n 33 

n oi 

n 02 

n 03 

n 12 
n 13 

n 23 



tt22 _ _ , ,2 

11 — £o^ p 



Am 2 Am 2 

,2 



U! 



+ 



n 2 



m 2m Am 2 



+ 



+ 



LU 2 -il 2 2m Am 2 2m{uj 2 - Vl 2 ) 

■ 2 i 2 7,2 



fc 2 



+ 7T 1 + 



exp 
fc 2 



exp 



2eB 



2eB 



1 H 1 H — H — + . 

2m 4m 2 2eB uj 2 - Q 2 \ m 2m Am 2 



exp 



2eB 



s u>. 



exp - 



2eB 



-e w. 



2m w 2 - Si 2 

fc||/c 2 



exp 



2eS 



p m(u> 2 — n 2 ) 



1 



= -te Lu 



fl 2 



U) 2 — fl 2 



exp 



2eB 



= -£oW 



, A; j. fcii 



p 



m 
fc±fc|| 
m 



1 

4m 
1 



2msl e 
1 



exp 



2eB 
T(w 2 + 2 v- 



+ 



2(w - n e ) (cj 2 - SI 2 ) 2 
1 2Tw!l e 



4m 2(w-fi e ) (w 2 -Sl 2 ) 2 



exp 
exp ( - 



A. 

2eB 
k\ 



2eB 



(BIO) 



14 



[1] G.G. Raffelt, Stars as Laboratories for Fundamental Physics (University of Chicago Press, Chicago, 1996). 

[2] J.B. Adams, M.A. Ruderman, and C-H. Woo, Phys. Rev. 129, 1383 (1963). 

[3] M.H. Zaidi, Nuovo Cimcnto 40, 502 (1965). 

[4] J.C. D'Olivo, J.F. Nieves, and P.B. Pal, Phys. Rev. D 40, 3679 (1989). 

[5] E. Braaten and D. Segel, Phys. Rev. D 48, 1478 (1993). 

[6] D.V.Gal'tsov and N.S.Nikitina, Zh. Eksp. Teor. Fiz. 62, 2008 (1972) [Sov. Phys. JETP 35, 1047 (1972)]. 

[7] V.V. Skobolev, Zh. Eksp. Teor. Fiz. 71, 1263 (1976) [Sov. Phys. JETP 44, 660 (1976)]. 

[8] L.L. DeRaad, Jr., K.A. Milton, and N.D. Hari Dass, Phys. Rev. D 14, 3326 (1976). 

[9] A.N. Ioannisian and G.G. Raffelt, Phys. Rev. D 55, 7038 (1997). 

[10] V. Canuto, C. Chiuderi, and C.K. Chou, Astophys. Space. Sci. 7, 407 (1970). 

[11] V. Canuto, C. Chiuderi, and C.K. Chou, Astrophys. Space. Sci 9, 453 (1970). 

[12] D.B. Melrose and A.J. Parle, Aust. J. Phys. 36, 755 (1983). 

[13] D.B. Melrose, Aust. J. Phys. 36, 775 (1983). 

[14] D.B. Melrose and A.J. Parle, Aust. J. Phys. 36, 799 (1983). 

[15] A. A. Sokolov and I.M. Ternov, Synchrotron Radiaton (Pergamon Oxford, 1968). 

[16] J.F. Nieves, in Proceedings of the IV Mexican School of Particles and Fields, edited by J. Lucio and A. Zepeda (World 

Scientific, Singapore, 1992), p. 283-319. 

[17] D.B. Melrose, Plasma Phys. 16, 845 (1974). 

[18] G.I. Svetozarova and V.N. Tsytovich, Izv. Vysshikh. Uchebn. Zaved. Radiofiz. 5, 658 (1962). 

[19] D.B. Melrose and R.J. Stoneham, J. Phys. A 10, 1211 (1977). 

[20] D.B. Melrose, Instabilities in Space and Laboratory Plasmas, (CUP, Cambridge 1986). 

[21] J.G. Kirk and N.F. Cramer, Aust. J. Phys. 38, 715 (1985). 

[22] T.H. Stix The Theory of Plasma Waves, (McGraw Hill: New York, 1962). 

[23] Y. Kohyama, N. Itoh, and H. Munakata, Astrophys. J. 310, 815 (1986). 

[24] G.G. Pavlov, Yu.A. Shibanov, and D.G. Yakovlev, Astrophys. Space. Sci. 73, 33 (1980). 

[25] P. Elmfors, D. Grasso, and G.G. Raffelt, Nucl. Phys. B 479, 3 (1996). 

[26] I.S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Series and Products, (Academic Press, 1980), 8.970. 



15 



FIG. 1. The W diagram, Z diagram and the V-A diagram for the plasma neutrino process 



FIG. 2. Comparison of emission from the ordinary mode at 8 = ir/2, both with and without the axial vector coupling. The 
power is in units of Wm~ 3 and the electron number density is in units of m -3 . The temperature is 10 s K and the magnetic 
field is 0.1 B c . 

FIG. 3. Temperature and magnetic field regimes for which the axial vector coupling is likely to be important in astrophysical 
objects. 
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TABLE I. The natural modes of a cold plasma 



Modes at 9 = 0, k = (0, 0, 1) 



Mode 



Dispersion Relation 



Polarization Vector 



Longitudinal 
Ordinary 

Extraordinary 



Modes at = tt/2, k = (1,0,0) 



Mode 



P = 



n 2 M = S + D 



n 2 M = S-D 



Dispersion Relation 



e= (0,0,1) 



e= 71 (M,0 ) 



e=^(l,-i,0) 



Polarization Vector 



Ordinary 



Extraordinary 



«M = P 

S 2 -D 2 



e= (0,0,1) 
_ (A-AS.0) 

yi) 2 + s 2 



17 



V 



W 



W Diagram 

v 



Z Diagram 

v 



4- Vertex Diagram 
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